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ABSTRACT 

The most usual tracer of molecular gas is line emission from CO. However, the re- 
liability of that tracer has long been questioned in environments different from the 
Milky Way. We study the relationship between H2 and CO abundances using a fully 
dynamical model of magnetized turbulence coupled to a chemical network simplified 
to follow only the dominant pathways for H 2 and CO formation and destruction, and 
including photodissociation using a six-ray approximation. We find that the abun- 
dance of H2 is primarily determined by the amount of time available for its formation, 
which is proportional to the product of the density and the metallicity, but insensitive 
to photodissociation. Photodissociation only becomes important at extinctions under 
a few tenths of a visual magnitude, in agreement with both observational and prior 
theoretical work. On the other hand, CO forms quickly, within a dynamical time, but 
its abundance depends primarily on photodissociation, with only a weak secondary 
dependence on H2 abundance. As a result, there is a sharp cutoff in CO abundance at 
mean visual extinctions Ay < 3. At lower values of Ay we find that the ratio of H2 col- 
umn density to CO emissivity Xqq, oc Ay 3 ' 5 . This explains the discrepancy observed 
in low metallicity systems between cloud masses derived from CO observations and 
other techniques such as infrared emission. Our work predicts that CO-bright clouds 
in low metallicity systems should be systematically larger or denser than Milky Way 
clouds, or both. Our results further explain the narrow range of observed molecular 
cloud column densities as a threshold effect, without requiring the assumption of virial 
equilibrium. 
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1 INTRODUCTION 

Observed star formation takes place within giant molecular 
clouds (GMCs), so understanding how these clouds form and 
evolve is a key step towards understanding star formation 
on galactic scales. The main chemical constituent of any 
GMC is molecular hydrogen (H2). However, it is extremely 
difficult to directly observe this molecular hydrogen in situ. 
Radiative transitions in H2 are weak, owing to the molecule's 
lack of a permanent dipole moment. Moreover, the lowest 
lying rotational energy levels of H2 are widely spaced, and so 
are very rarely excited in gas with the temperatures typical 
of GMCs, T ~ 10-20 K. For this reason it is common to 
use emission from carbon monoxide (CO), the second most 
abundant molecule in GMCs, as a proxy for the H2. 

In order to use CO as a proxy for H2, however, it is 
necessary to understand the relationship between the dis- 
tributions of these two molecules. Both are readily dissoci- 



ated by the absorption of ultraviolet photons with energies 
below the Lyman limit of atomic hydrogen, so in low den- 
sity, low extinction gas such as the warm neutral component 
of the interstellar medium (ISM), the abundances of both 
molecules will be small. However, the two molecules form 
rather differen tly: H2 forms predomin antly on the surface 
of dust grains l|Gould fc Salpeterll 1963th while CO forms al- 
most exclusively in the gas phase, via any one of a number 
of chains of ion-neutral or n eutral-neutral reactions (see e.g. 
ISternberg fc Dalgarndll995l for a useful summary of CO for- 
mation chemistry). Moreover, H2 can protect itself from ul- 
traviolet radiation via self-shielding, whic h becomes effective 
for r elatively low H2 column densities l|Draine k, Bertoldl 
1 1996T) . The corr esponding process for CO is less effective 
( Lee et alj|l996h . requiring a higher column density, and 
the low abundance of carbon relative to hydrogen in the 
ISM means that the required column density typically cor- 
responds to a situation in which a large fraction of the avail- 
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able carbon is already locked up in the form of CO. For these 
reasons, we would expect the CO/H2 ratio to vary through a 
cloud, and this expectation is conf irmed by detailed models 
of slab-like or spherical clouds (e.g. lRollig et al.ll2007l ). These 
find that the transition from atomic hydrogen to molecular 
hydrogen occurs at a point closer to the cloud surface than 
the transition from ionized carbon, via neutral atomic car- 
bon, to CO. 

Although spatial variations in the CO/H2 ratio within 
a CMC would appear to make it difficult to use CO as a 
proxy for H2, observations of Galactic GMCs show that in 
fact there appears to be a good correlation between the inte- 
grated intensity of the J = 1 — > rotationa l transition line 
of 12 CO and the H2 column density (see e.g. Dickman 19781; 
Sanders et all 1 1984 ISolomon et all 1 19871 ; Istrong fc Mattoxl 



19961 ; iDame et al.l 120011 ). A number of independent studies 



have shown that GMCs in the Galactic disk all have CO-to- 
H2 conversion factors that are approximately 
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where Wco is the velocity-integrated intensity of the CO 
J = 1 — > emission line, averaged over the projected area 
of the GMC, and Nn 2 is the mean H2 column density of 
the GMC, averaged over the same area. Although the for- 
mer can be directly observed, the latter cannot. However, 
for nearby clouds, it can be inferred from measurements of 
the diffuse 7-ray flux produced by interactions between high 
energy cosmic rays and atomic hydrogen, atomic helium and 
H2. The 7-ray flux along a given line of sight depends on the 
total hydrogen column density along that line of sight. Since 
the atomic hydrogen column density can be measured via its 
21 cm emission, the H2 column density can be inferred. As 
the values for Xco obtained in this way are consistent with 
those obtained by assuming that the GMCs are in virial 
equilibrium and using the obser ved linewidth-size re lation 
to compute the cloud mass (e.g. ISolomon et ail 1 19871 ). it is 
generally accepted that CO emission is indeed a good proxy 
for H2 mass in nearby GMCs. 

However, the issue of the environmental dependence of 
Aco remains highly contentious. In gas with a lower metal- 
licity, or a higher ambient UV radiation field, CO pho- 
todissociation will be more effective, and the amount of 
CO in the cloud will be smaller, with the CO-rich gas oc- 
cupying a smaller volu me than in the Galactic case (e.g. 
iMalonev fc Blacklll988l . Molina et al., in prep.). The mean 
H2 abundance may also be smaller, but the greater role 
played by H2 self-shielding means that we do not expect 
the H2 abundance to be nearly as sensitive to changes in 
the metallicity. Because of this, there are good theoretical 
reasons to expect the relationship between CO emission and 
H2 mass to change as we change the metallicity. However, 
observational efforts to test this yield inconsistent results. 

Measurements of Ago that assume that extragalac- 
tic GMCs are in virial equilibrium and use a virial anal- 
ysis to determine the cloud mass generally find values for 
Ago that are similar to those obtained in th e Milky Way , 
with at most a weak metallicity dependenc e (j Wilsonl 1 19951 ; 
iRosolowskv et all 1200 j ; iBolatto et ail 120081 ). On the other 
hand, measurements that constrain GMC masses using other 
techniques that do not depend on the CO emission, such 
as by measuring the far-infrared dust emission, consistently 



find values for Ago that are much larger than the Galac- 
tic value a nd that are suggestive of a strong metallicity de- 
pend ence <|lsraellll997l ; iRubio et al. 1 12004 iLerov et al.ll2007l . 

Numerical simulations provide us with one way to ad- 
dress this observational dichotomy. If we can understand 
the distribution of CO and H2 in realistic models of GMCs, 
then we may begin to understand why the different types 
of observation give such different results. However, until 
very recently, the ability of simulations to address this is- 
sue has been quite limited. Sophisticated treatments of gas 
and grain chemistry, radiative heating and cooling, and ra- 
diative transfer have bee n developed to m odel GMCs (see, 
for instance, the paper bv lRollig et ai1l2007l , which compares 
results from a number of popular codes), but the very com- 
plexity of these models limits their applicability to highly 
simplified geometries: typically the adopted cloud models 
are one- dimensional, assuming either s pherical symmetry 
(e.g. IStoerzer. Stutzki fc Sternberg||l996l ) or a semi-infinite, 
uniform slab (e.g. iLe Petit et al.ll2006h . 

Moreover, the modelling often assumes that the chem- 
istry of the clouds is in equilibrium, which is valid only if 
all of the chemical timescales are much shorter than any 
dynamical timescale associated with the GMCs. Since real 
GMCs are observed to be highly inhomogeneous, and to 
be dominated by super sonic turbulent motions (see e.g. 
iMac Low fc Klessenl |200"3 . and references therein), the ap- 
plicability of the results from these simple spherical or slab 
models to real clouds is open to question. On the other 
hand, previous attempts to accurately model the turbu- 
lent dynamics of GMCs, and the inhomogeneous and in- 
termittent density structure that is created by this tur- 
bulence typically have avoided modelling the cloud chem- 
istry, rendering them also of limited use in addressing this 
question. The few studies that have attempted to model 
bot h the turbulence and the chemistry self-consistently 
Joulain et al.lll998l; iLesaffre. Gerin fc Hennebelldl2007l; 



iGodard et al.l 20091) have typically done so b y reducing the 
dimen sionality of the problem. For instance, IGodard et alj 
(2009) present results on the chemical evolution of gas 
passing through magnetized, turbulent vortices, using two- 
dimensional simulations of these vortices, and then construct 
models for lines of sight through diffuse clouds by summing 
up the contributions from a number of different vortices. 
However, it is unclear whether the results obtained in this 
fashion are the same as those that would be obtained using 
a fully three-dimensional model for the turbulence. 

In a previous paper (|Glover et alj|201ol . hereafter Paper 
I), we presented the first results from a project that aims to 
combine detailed chemical modelling with three-dimensional 
magnetohydrodynamical turbulence simulations in order to 
self-consistently model both the chemistry and the turbu- 
lent dynamics of the gas within a GMC. We showed that 
it is now computationally feasible to attempt this. Provided 
one makes a few simplifying assumptions regarding the ex- 
tent of the chemistry to be treated, and the treatment of 
the UV radiation field, it is possible to model both the H2 
and the CO chemistry of a GMC with acceptable accuracy 
within a moderate resolution dynamical simulation. Paper I 
presented results from a few trial simulations performed us- 
ing only one set of cloud properties (mean density, metallic- 
ity, etc). In this paper, we present results from a much larger 
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set of simulations that examine the sensitivity of H2 and CO 
formation to changes in the mean densities and metallicities 
of GMCs, and explore the consequences of this for the CO- 
to-H2 conversion factor and the dynamical structure of the 
observed clouds. 

In section [21 we discuss our numerical approach and the 
initial conditions used for our simulations. In section O we 
present the main results of our simulations, and use them 
to derive the approximate dependence of the CO-to-H2 con- 
version factor on cloud properties. In section [4j we discuss 
the major consequences of our findings. In section[5l we dis- 
cuss a few potential caveats regarding our approach, and 
show why they are unlikely to significantly affect our main 
results. Finally, in section[6l we close with a brief summary. 



2 SIMULATIONS 
2.1 Numerical method 

Our simulations were performed using a mod ified version of 
the ZEUS-MP ma gnetohydrodynamical code (|Normanl2000l ; 
lHaves et al.ll2006l ) fully described in Paper I. Our modifica- 
tions include the addition of a simplified treatment of hy- 
drogen, carbon and oxygen chemistry, composed of 218 reac- 
tions between 32 chemical species, together with a detailed 
atomic and molecular cooling function. This treatment in- 
cludes the formation of H2 on the s urface of dust grains, 
follow ing a prescription taken from iHollenbach fe McKed 
(1979). The formation rate of H2 per unit volume due to 
grain surface reactions is written as 



Rh 2 = 7h 2 



nnu cm 



(2) 



where n is the number density of hydrogen nuclei, riH is the 
number density of atomic hydrogen, Z is the metallicity of 
the gas, and 7h 2 is the rate coefficient for H2 formation, 
given by 



7h 2 
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where T2 = T/100K is the gas temperature in units of 100 K, 
Td,2 = Td/100 K is the dust temperature in units of 100 K, 
and where 



f: 



(4) 



l + 10 4 exp(-600/T d )' 

For simplicity, we fix the dust temperature at 10 K in all 
of the simulations presented in this paper. Aside from this 
reaction, we include no other grain surface chemistry in our 
current models. 

We model the photodissociation of H2 using prescrip- 
tions for dust shi e lding and H2 self- shielding taken from 
iDraine fe BertoldH |l996). For CO photodissociation, we ac- 
count for dust shielding, CO self-shielding and shielding by 
H2, using shielding functions taken from iLee et alJ (| 19961 ) 
for the last two terms. To compute the column densities of 
dust, H2 and CO required by this treatment of photodis- 
sociation, we use the six-r ay approximation introduced in 
Idover fc Mac Low! (|2007al ) . In this approximation, we com- 
pute photochemical rates in each zone in our simulation vol- 
ume by averaging over the six lines of sight that lie parallel 
to the coordinate axes. 



2.2 Initial conditions 

All of our simulations begin with an initially uniform gas 
distribution, with mean hydrogen nucleus number density 
no- The initial temperature of the gas is 60 K. However, be- 
cause the cooling time of the gas is significantly shorter than 
the dynamical time in all but our lowest density runs, the 
gas temperature quickly adjusts itself until the gas is close 
to thermal equilibrium, and hence our results are largely 
insensitive to our choice for the initial temperature. The ini- 
tial velocity field is turbulent, with uniform power between 
wavenumbers k = 1 and k = 2, and with an initial rms veloc- 
ity of 5 km s _1 . We drive the turbulence so as to maintain 
approximately the same rms velocity thro ughout the sim- 
ulatio n, fol lowing the meth od described in iMac Low et al.l 
l|l998l ) and lMac Lowl l|l999l ). 

We assume that the gas is magnetized, with an initially 
uniform magnetic field strength Bo = 5.85 piG and a field 
that is ini tially oriented paral l el to th e z-axis of the simu- 
lation. In I Glover fc Mac Low! (|2007bl ) we showed that the 
timescale for H2 formation in a turbulent, magnetized cloud 
is relatively insensitive to the value chosen for the magnetic 
field strength, and hence we do not expect our results to 
be particularly sensitive to our choice of Bo- The particu- 
lar value that we have chosen is the same as we adopted 
in Paper I, and is motiv ated by the 21 cm observations by 
iHeiles fc Trolandl l|2005l ) that find a median magnetic field 
strength of ~ 6/xG for the cold neutral medium. 

We assume that the gas has a uniform metallicity Z, and 
that the initial fractional abundances (by number) of carbon 
and oxygen relative to hydrogen are given by x c + = a;c,tot 
I . ', I x'lO- 4 (Z/Z ) and x = xo,tot = 3.16 x 10- 4 (Z/Z Q ), 
where a;c,tot and a;o,tot refer to carbon and oxygen in all 
forms (ionized, neutral, or incorporated into molecules). For 
solar metallicity, this corresp onds to the values given i n 
ISembach efaZI (|2000i ) . As in Idover fc Mac Low! (|2007tJ ). 
we ignore the effects of metals other than C or O, in order 
to minimize the computational requirements of our simula- 
tions. We assume a standard dust-to-gas ratio for our solar 
metallicity runs, and also assume that in our lower metallic- 
ity runs, the dust-to-gas ratio scales linearly with metallicity. 
We assume that the dust has an extinction curve character- 
ized by i?v = 3.1, and that this is independent of metallicity. 
We also fix the dust temperature at 10 K in every run. 

We adopt a fixed cosmic ray ionization rate £ = 
10~ 17 s _1 . For the majority of our simulations, we use the 
same spectral shap e and normal ization for the interstellar 
radiation field as in lDrainel (|l978l ). For brevity, when we re- 
fer to the strength of the radiation field in a simulation, we 
typically give it in terms of Go, the ratio of the field strength 
relative to the strength of the Draine field. Thus, in the ma- 
jority of our simulations, Go = 1.0. We also performed a 
few simulations without a radiation background, i.e. with 
Go = 0.0. The effects of increasing Go to much larger values 
are examined in a separate paper (Glover et al., in prep). 

We perform our simulations in a cubical volume of side 
L, to which we apply periodic boundary conditions. As we 
have previously noted in Paper I, the use of periodic bound- 
ary conditions for the gas together with what are effectively 
non-periodic boundary conditions for the radiation is not 
self-consistent, and we adopt this arrangement purely for 
the computational convenience that it offers. We do not in- 
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Table 1. List of simulations 



ID 


no (cm °) 


Z/Z Q 


Box size (pc) 


G 


n30 


30 


1.0 


20 


1.0 


nlOO 


100 


1.0 


20 


1.0 


nl80 


180 


1.0 


20 


1.0 


n300 


300 


1.0 


20 


1.0 


nlOOO 


1000 


1.0 


20 


1.0 


n300-Z06 


300 


0.6 


20 


1.0 


nl00-Z03 


100 


0.3 


20 


1.0 


n300-Z03 


300 


0.3 


20 


1.0 


n30-Z01 


30 


0.1 


20 


1.0 


nlOO-ZOl 


100 


0.1 


20 


1.0 


n300-Z01 


300 


0.1 


20 


1.0 


nlOOO-ZOl 


1000 


0.1 


20 


1.0 


nl000-Z003 


1000 


0.03 


20 


1.0 


n30-L5 


30 


1.0 


5 


1.0 


nlOO-L5 


100 


1.0 


5 


1.0 


n300-L5 


300 


1.0 


5 


1.0 


nl000-L5 


1000 


1.0 


5 


1.0 


n30-UV0 


30 


1.0 


20 


0.0 


nlOO-UVO 


100 


1.0 


20 


0.0 


n300-UV0 


300 


1.0 


20 


0.0 


nlOOO-UVO 


1000 


1.0 


20 


0.0 



Note: Go is the ultraviolet field strength in units of the iDraind 
l|l978l l field 

elude the effects of self-gravity in our current simulations, 
but intend to explore them in future work. 

We have performed simulations with a range of different 
initial densities no and metallicities Z, as detailed in Tabled] 
For most of these simulations, we set L = 20 pc, but we have 
also performed a few simulations with L = 5 pc. For all of 
these simulations, we adopt a numerical resolution of 128 3 
zones, which we showed in Paper I was sufficient to yield 
converged results for the mass-weighted mean abundances 
of CO and H2, at least to the level of accuracy required 
here. 

We evolved all of our runs until f cnc i = 1.8 x 10 14 s ~ 
5.7Myr, corresponding to roughly three eddy turnover times 
for the turbulence. In Paper I, we showed that this gave 
sufficient time for the density distribution, temperature dis- 
tribution and the chemical abundance of CO all to reach 
a statistical steady state. We expect the same to be true 
for the majority of runs examined here. However, as we will 
see later, there is evidence that the CO abundance in some 
of our lowest density runs has not reached a steady state 
by t = tend- Therefore, in four of our lowest density runs 
(runs n30, nlOO, n30-Z01 and nlOO-ZOl), we did not stop 
at this point, but continued to evolve them until a time 
text = 7.2 x 10 14 s ~ 22.8 Myr. We note, however, that it 
is questionable whether real molecular clouds could live for 
so long without undergoing localized gravitational collapse 
and forming stars, physics that is not included in our current 
models. 



3 RESULTS 

3.1 Density and temperature structure of the gas 

Before considering the molecular abundances produced in 
the various simulations, it is perhaps useful to briefly dis- 



cuss the density and temperature distributions generated 
by the turbulence. As an illustrative example, we plot in 
Figure [IJl the mass-weighted density probability distribu- 
tion function (PDF) for runs n300 and n30. The PDFs are 
plotted in terms of the dimensionless logarithmic density 
contrast s = ln(p/po), where po is the mean density of the 
gas. It is clear from the plot that the two PDFs have a simi- 
lar form when plotted in dimensionless units. Both PDFs are 
roughly log-normal, with similar widths and means, as ex- 
pected from p revious studies of su personic turbulence in in- 
terstellar gas (|Padoan et aZ.|[l"997t ). Deviations from the log- 
normal form are apparent in the wings of the PDFs, where 
one would expect t he effects of turbulent intermittency to be 
most pronounced (TKritsuk et al.ll2007l ; ISchmidt et al.ll2009l : 
iFederrath et al.ll2010r i.~ 

It is also interesting to compare the temperature PDFs 
of the same two runs (Figure lb). If we do so, we find that 
the gas is systematically warmer in run n30 than in run n300, 
owing to the increased effectiveness of photoelectric heating 
and reduced effectiveness of CO cooling in the lower density 
run. We find similar behaviour in our other runs, and so for 
reasons of brevity, we will not plot their density and temper- 
ature PDFs here. In each case, we find a roughly log-normal 
density PDF, as expected, with a width that is correlated 
with the mean temperature of the gas. Reducing the den- 
sity and/or the metallicity of the gas tends to increase its 
mean temperature, both by reducing the efficiency of cool- 
ing and by increasing the influence of photoelectric heating, 
and hence the lower density and lower metallicity runs have 
narrower density PDFs than the higher density or higher 
metallicity runs. This trend is also apparent if one looks at 
the clumping factor of the gas, defined as 



C = 



(£!> 
(P) 2 ' 



(5) 



which ranges from C ~ 4.5 in run nl000-Z003 to C ~ 6.8 
in run nlOOO, with the values from the other runs spread 
fairly evenly between these two extremes. These values are 
consistent with those inferred fr om models of some Galac- 
tic photodissociation regions (e.g. lKramer et al1l2004l . whose 
results imply a mean value for C of roughly 5 in W3 Main) , 
although other regions show evidence for much higher values 
of C (e.g. iFalgarone fc Phillipslll996r ). 

Although these resul ts are not particu l arly su rpris- 
ing given our findings in Idover fc Mac Low! (l2007bh and 
in paper I, they do serve to demonstrate that the den- 
sity structure produced by the turbulence is rather dif- 
ferent from the structure that is often assumed in mod- 
els of clumpy photodissociation regions (PDR). These 
models commonly assume that the density structure of 
a PDR can be represented by a population of clumps 
(usually assumed to be spherical, but with densities and 
sizes that can vary from clump to clump), embedded 
in some much lower density interchrmp medium (see 



e.g. iBurton, Hollenbach & Tielensl 


1990;[Meixner & Tielensl 


1993; Spaans & van Dishoeck! 1997 


; Mookeriea et all 2006). 



However, the density distributions generated by the turbu- 
lence in our simulations do not fit comfortably into this pic- 
ture. The fact that we find density PDFs that are smooth, 
single-peaked functions without even a hint of bimodal be- 
haviour suggests that any attempt to partition the gas into 
high density clumps embedded in a much lower density 



H 2 Versus CO Abundance 5 




Table 2. Mass-weighted mean abundances of H2 and CO at our 
default stopping time, t = t cn( j 
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Figure 1. (a) Mass- weighted PDF of the logarithmic density 
contrast s = ln(p/po)i where po is the mean density of the gas, 
plotted for runs n300 (solid line) and n30 (dashed line). Both 
of the PDFs display the log-normal shape characteristic of quasi- 
isothermal supersonic turbulence, (b) Mass- weighted temperature 
PDF for runs n300 (solid line) and n30 (dashed line). The char- 
acteristic temperature of the gas in run n30 is much higher than 
in run n300, owing to a combination of the effects of the lower 
gas density and the low CO abundance (see A3. 21 below), which 
render cooling much less effective, and the lower level of dust 
extinction, which significantly increases the importance of pho- 
toelectric heating. Note that while the relatively sharp cutoff in 
the PDF for run n30 at T < 30 K appears to be real, the corre- 
sponding feature at T < 10 K in the PDF for run n300 is not; the 
latter is a consequence of the fact that for technical reasons, our 
cooling function cuts off at T = 10 K (see Paper I). 



medium will inevitably wind up being somewhat arbitrary, 
since there is no particular density or range of densities that 
one can pick out of the PDF as corresponding to a 'clump', 
and hence no obvious way to determine where a 'clump' ends 
and the interclump gas begins. This unfortunately makes it 
difficult to compare our results with those coming from these 
clumpy PDR models, but does also serve to emphasize the 
great importance of modelling the dynamics and chemistry 
of the gas in a self-consistent fashion. 



ID 


(zh 2 )m 


(zco)m 




n30 


0.351 


9.31 


-9) 




nlOO 


0.689 


2.97 


-6) 




nl80 


0.847 


1.92 


-5) 




n300 


0.943 


6.41 


-5) 




nlOOO 


0.998 


1.27 


-4) 






0.874 


2.10 


-5) 




nl00-Z03 


0.357 


2.19 


-8) 




n300-Z03 


0.712 


2.74 


-6) 




n30-Z01 


0.027 


2.27 


-12) 




nlOO-ZOl 


0.136 


8.27 


-11) 




n300-Z01 


0.385 


1.66 


-7) 




nlOOO-ZOl 


0.816 


2.32 


-6) 




nl000-Z003 


0.449 


6.08 


-8) 




n30-L5 


0.342 


3.80 


-10) 




nl00-L5 


0.817 


1.26 


-8) 




n300-L5 


0.982 


9.00 


-6) 




nl000-L5 


0.998 


8.17 


-5) 




n30-UV0 


0.345 


1.84 


-5) 




nlOO-UVO 


0.687 


8.39 


-5) 




n300-UV0 


0.944 


1.39 


-4) 




nlOOO-UVO 


0.998 


1.41 


-4) 




Table 3. Mass 


-weighted 


mean 


abundances of H2 


end of our extended runs 








ID (zh 2 )m at t cx t 


(^co)m 


&t text 


n30 


0.774 




2.63( 


-8) 


nlOO 


0.972 




1.55( 


-5) 


n30-Z01 


0.055 




6.78(- 


-12) 


nlOO-ZOl 


0.445 




2.84(- 


-10) 



3.2 Molecular abundances: dependence on 
extinction and gas density 

Let us now turn our attention from the density and temper- 
ature of the gas to its chemical content. In Table [2] we list 
the mass- weighted mean abundances of H2 and CO, denoted 
(ih 2 )m and (:tco}m respectively, at our default end-point 
t = tend = 5.7 Myr, for each of the simulations that we have 
run. These mean abundances are defined as 

/ \ 2 AU V n H 2 («, j, k) 



(6) 



and 



<*CO>M=^:Ay£ 



nco(i,j, k) 
n{i,j, k) 



p(i,j,k), 



(7) 



respectively, where M is the total mass in the simulation, 
AV is the volume of a single grid zone, n(i,j, k) is the num- 
ber density of hydrogen nuclei in zone (i,j,k), nu_ 2 (i, j,k) 
and nco(i,j, k) are the H2 and CO number densities in the 
same zone, p(i,j,k) is the mass density in that zone, and 
where we sum over all grid zones. Note that the inclusion 
of the factor of two in our definition of (ih 2 )m is to en- 
sure that when the hydrogen is fully molecular, we will have 
(£h 2 }m = 1.0, rather than 0.5, as this will hopefully be less 
confusing for the reader. 
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Figure 2. (a) Mass- weighted mean H2 fraction, plotted as a func- 
tion of the mean visual extinction (Ay). Different symbols are 
used for runs with different metallicities, box sizes or UV field 
strengths, as outlined in the legend. The runs performed using 
an L = 5 pc box (blue circles) produce as much H2 as the corre- 
sponding L = 20 pc runs (crosses), and so the former are simply 
displaced to the left of the latter by a factor of four. This fact, 
together with the fact that the runs performed without a UV 
background (red stars) follow the same apparent correlation as 
the runs with a UV background, provide convincing evidence that 
the apparent correlation between (ih 2 )m and (Ay) is artificial: 
the true correlation is between (:eh 2 )m an< ^ n oZ, and the appar- 
ent correlation with (Ay) arises because (Ay) oc n^LL. (b) As 
(a), but for the mean mass- weighted CO fraction. In this case, the 
L = 5 pc runs show the same correlation between the mean CO 
fraction and the mean extinction as the L = 20 pc runs, but the 
runs performed without a UV background do not. This suggests 
that the correlation between (ico)m an d (Ay) is real. 



In Figured, we plot the values of (ih 2 )m at t = t cn d 
as a function of the mean visual extinction, (Av), and in 
Figure [2p, we give a similar plot of the values of (xco)m at 
t = t cn d- The mean visual extinction (Av) is defined as 



(8) 



3,k 



where TV is the number of grid zones along one edge of the 
box, Ay.j.fc is the visual extinction along the line of sight 
running parallel to the a;-axis and with y and z coordinates 
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Figure 3. (a) Mass-weighted mean H2 fraction, plotted as a 
function of mean density no . We see a clear dependence on den- 
sity at fixed metallicity, and a clear dependence on metallicity at 
fixed density, consistent with (xjj 2 )m correlating with the prod- 
uct of density and metallicity. (b) As (a), but for the mean mass- 
weighted CO fraction. Again, the CO fraction depends on density 
at fixed metallicity, and depends on metallicity at fixed density, 
but now the CO fraction also depends strongly on the box size L. 
Density by itself, without knowledge of L or Z, predicts the CO 
fraction poorly. 



(j, k), and where we sum over all sets of (j, k) corresponding 
to zone centres. Note that our choice to project along the 
i-axis is arbitrary, and that projecting along the y or z axes 
would yield the same value of (Av), which indeed is simply 
given by 



(Av) = fconvUoL, 



(9) 



where / C onv is the conversion factor between the column den- 
sity of H nuclei and the visual extinction. This is given by 
/conv = 5.348 x 10~ 22 (Z/Z Q ) mag cm 2 for our choice of dust 
properties. 

To understand these results, first consider only the runs 
performed with our default values for L and Go- Figure [2^ 
shows that if we disregard the L — 5 pc runs, then a clear 
correlation exists between (a;H 2 )M and (Av), with larger val- 
ues of (Av) yielding larger values of (a;H 2 }M- Similarly, if we 
disregard the Go = 0.0 runs, then Figure [2p shows that 
there is also a clear, and far stronger, correlation between 
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(xco)m and (Av). If we compare the runs with (Ay) = 1 
and (Ay) — 10, then we see a factor of two to three change in 
(:eh 2 }m, compared with a change in (ico)m of three to four 
orders of magnitude. Plotting (s;h 2 }m and {xco)m against 
the mean gas number density no shows no clear correla- 
tion in either case (Figure [3]), if we consider the data as a 
whole, but also shows that both (:eh 2 }m and {xco)m corre- 
late well with density at fixed metallicity, and with metal- 
licity at fixed density. 

At first glance, these results suggest that both the H2 
and the CO abundances are controlled by photodissocia- 
tion, with the CO responding far more strongly to changes 
in (Ay) than the H2. However, the results of the L — 5 pc 
runs challenge this interpretation. They lie where we would 
expect in the plot of (ico)m versus (Ay) (Figure Hp), but 
they lie to the left of the corresponding L — 20 pc runs in 
the plot of (ih 2 )m versus (Ay) (Figure [2^,). Indeed, shift- 
ing them uniformly to the right by a factor of four brings 
them into good agreement with the other runs. This result 
can be understood if we realize that the real correlation re- 
sponsible for the behaviour seen in Figure [2K is riot one be- 
tween (:eh 2 }m and (Ay), but is instead a correlation between 
(ih 2 )m and noZ. Since (Ay) oc noZ, we find an apparent cor- 
relation between (xh 2 }m and (^4v) if we compare only runs 
with the same box size L. However, if we change L without 
changing noZ, we will still obtain the same value of (:th 2 }m, 
but find that it now corresponds to a different mean visual 
extinction. 

This result strongly suggests that the differences in the 
mean H2 abundance found at the end of the different runs 
are not a result of photodissociation, and that this plays 
little or no role in regulating the amount of H2 present in 
the gas. Instead, the dependence of (ih 2 )m on noZ suggests 
that the amount of H2 present depends simply on the H2 
formation time, tform,H 2 , which scales as tf or m,H 2 oc l/(noZ). 
If noZ is small, then the H2 formation time is long, and the 
hydrogen in the gas simply does not have sufficient time to 
become fully molecular by the end of the simulation, de- 
spite the boost to the H2 formation rate resulting from the 
presence of the turbulent compressions. On the other hand, 
when noZ is large, the hydrogen can easily become almost 
fully molecular. 

In contrast, the CO abundances in these runs respond 
in the fashion that we would expect them to if they were pri- 
marily controlled by photodissociation. Decreasing L while 
holding uqZ constant leads to a sharp drop in the CO abun- 
dance, implying that in this case the correlation with (Ay) 
is real, and not simply a consequence of the relationship 
between noZ and (Ay). 

We can confirm this interpretation by considering what 
happens in the runs that are performed with Go = 0.0 to 
eliminate all photodissociation. We see from Figure [2ji that 
in this case, we still recover the same correlation between 
{xh 2 )m and (Ay). However, in the absence of an ultraviolet 
background, the H2 abundance cannot depend on the vi- 
sual extinction, as the only physical processes in the simula- 
tion that depend on the visual extinction are those involving 
the attenuation of the ultraviolet background. Therefore, the 
correlation between (xh 2 )m and (Ay) cannot be real, and 
must simply be a consequence of the actual correlation of 
between (xh 2 )m and noZ, or rather the inverse correlation 
between (xh 2 )m and the mean H2 formation time. On the 



other hand, Figure [2p shows that the CO abundance at low 
(Ay) is much larger in the runs without a UV background, 
consistent with a picture in which the CO abundance is pri- 
marily controlled by photodissociation. 

Figure [2p also demonstrates that the CO abundance 
is not determined solely by photodissociation, since even 
in the absence of the UV background, we see a systematic 
decrease in the CO abundance as we move to lower densi- 
ties. This decrease is driven primarily by the decrease in the 
mean H2 abundance, since H2 is a key component in most 
of the chemical pathways capable of forming gas-phase CO. 
We have verified this by running simulations similar to runs 
n30-UV0 and nlOO-UVO that start with all of the hydrogen 
in molecular form. The CO abundances produced in these 
runs are significantly higher than in the runs that start with 
atomic hydrogen: in the lower density run, roughly 60% of 
the available carbon is converted to CO, while in the higher 
density run, roughly 95% is converted. In comparison, the 
fraction of carbon converted to CO in runs n30-UV0 and 
nlOO-UVO is 13% and 60%, respectively. Nevertheless, it is 
clear from Figure [5p that the dependence of the CO abun- 
dance on the H2 abundance is a relatively small effect in 
comparison with the dependence of the CO abundance on 
the mean visual extinction in the runs with a UV back- 
ground: the former is responsible for changes in (ico)m of 
a factor of a few, while the latter is responsible for changes 
of several orders of magnitude. 

We have thus arrived at the first key result of this pa- 
per: in a turbulent molecular cloud, the mean H2 abundance 
depends primarily on the time taken to form the H2 and is 
insensitive to photodissociation, while the CO abundance 
is determined primarily by photodissociation, with only a 
secondary dependence on the H2 fraction. 

3.3 Molecular abundances: time dependence 

We can get further confirmation of the results of the previous 
section by considering the time dependence of {xh 2 )m and 
{xco)m- In Figuref!] we plot the evolution of these quantities 
as a function of time for a set of runs with different mean 
densities at solar metallicity, while in Figure O we show a 
similar plot for runs with a metallicity Z = 0.1 Z@. Looking 
at the evolution of H2 fraction with time in these figures, we 
see that the time required to convert a large fraction of the 
initial atomic hydrogen to molecular hydrogen decreases as 
we increase the density or the metallicity, just as we would 
expect given the results of the previous section. Moreover, 
it is clear that in most of the runs, the mean H2 fraction 
has yet to reach equilibrium at the end of the run. Finally, 
extrapolation of the curves beyond six Myr suggests that the 
H2 fraction in all of the runs should eventually reach a value 
of order unity; there is no indication that the equilibrium 
value of (s;h 2 )m will be significantly smaller than one in any 
of the runs. 

The time evolution of (ico)m in the different runs 
is qualitatively different from that of (sh 2 )m- CO forms 
rapidly at the beginning of the simulations, but after roughly 
a million years, the rate of increase of the CO abundance 
slows or stops. The amount of CO formed during this ini- 
tial period of rapid evolution depends strongly on the mean 
visual extinction of the gas, with the runs with low (Ay) 
forming considerably less CO than the runs with high (^4v). 
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The subsequent evolution of the CO abundance depends on 
the value of (xh 2 }m at this point. If (zi^m ~ 1, as in runs 
n300 or nlOOO, then there is little further evolution of the 
CO abundance. On the other hand, if the mean H2 abun- 
dance remains significantly smaller than unity at this early 
point, then the CO abundance continues to increase, driven 
by the increasing availability of H2, and may double or triple 
by the end of the run. 

It is also evident that (xco)m displays far larger fluc- 
tuations in runs with low mean visual extinctions (e.g. runs 
n30 or nlOO-ZOl) than in runs with high mean visual extinc- 
tions. This is a consequence of the difference in the spatial 
distribution of the CO in a low extinction run compared to 
a high extinction run. When (Av) is small, the CO fraction 
in most of the gas is very small, owing to the lack of effec- 
tive shielding of the CO against UV photodissociation. In 
this case, most of the CO in the simulation is found in a 
small number of dense, comparatively well-shielded clumps. 
Since our simulations were performed without self-gravity, 
these clumps are transient objects, and are continually being 
disrupted and reformed by the turbulence. This continual 
variation in the number and nature of the CO-rich clumps 
leads to the observed variation in the mass-weighted mean 
CO abundance. On the other hand, when (Av) is large, far 
more of the gas has a high CO fraction, and many more 
clumps and other dense structures such as sheets and fila- 
ments contribute to (ico)m. Although all of these structures 
are transient, just as in the low (Ay) runs, the fact that a 
much larger number of structures contribute to {xco)m than 
in the low (Av) runs means that the variation in {xco)m 
resulting from the creation or destruction of any particular 
clump or filament is much smaller. 

In Figure[6l we show how {xh 2 )m and (xco)m vary with 
time in runs n30-UV0, nlOO-UVO, n300-UV0, and nlOOO- 
UVO, all performed without an ultraviolet background. The 
evolution of the H2 fraction in these runs is very similar to 
that in the runs with a UV background, again showing that 
photodissociation plays little or no role in regulating the 
growth of the H2 abundance. On the other hand, the evolu- 
tion of the CO abundance in these runs is clearly different 
from that in the runs with a UV background. The increase 
in the CO abundance in these runs more clearly tracks the 
increase in the H2 abundance, and there is no longer any sig- 
nificant feature at t ~ 1 Myr (other than in run nIOOO-UVO, 
which has converted almost all of its carbon into CO by this 
point). The CO abundances produced in the lower density 
runs are also clearly much larger than in the corresponding 
runs with a UV background, and there is no indication that 
they have yet reached equilibrium. 

Finally, in Figure we plot the time evolution of 
(ih 2 )m and {xco)m in the four runs that we evolved for a 
much longer time. This plot further emphasizes that time is 
the primary factor controlling the H2 abundance in the runs 
in which the product wqZ is small. For instance, at t ~ 5Myr, 
the mass- weighted mean H2 fraction in run nlOO-ZOl is only 
10%, but it has increased to roughly 40% by t ~ 20 Myr, i.e. 
it has increased almost linearly with time. The only one of 
these runs in which H2 photodissociation appears to play an 
important role is run n30-Z01, in which the mean H2 fraction 
appears to saturate at roughly 5%. The mean visual extinc- 
tion of the gas in this run is only 0.1 magnitudes, while in run 
nlOO-ZOl, which is not significantly affected by the UV, it is 
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Figure 4. (a) Evolution of the mass-weighted mean H2 abun- 
dance with time in runs nlOOO (solid line), n300 (dashed line), 
nl80 (dash-dot-dot-dotted line), nlOO (dash-dotted line) and n30 
(dotted line), (b) As (a), but for the mass- weighted mean CO 
abundance. The large fluctuations seen in run n30 occur because 
the dominant contribution to {xqo)m in this simulation comes 
from only a few transient CO-rich clumps. As the turbulence 
forms and destroys them, their contribution to (ico)m under- 
goes large fluctuations. In the higher density runs, many more 
structures contribute significantly to (icq)m an d so the influ- 
ence of the formation or destruction of any one clump is much 
smaller. 



approximately 0.3 magnitudes, suggesting that the critical 
value below which UV photodissociation begins to signifi- 
cantly limit the H2 abundance lies somewhere in the range 
(Av) ~ 0.1-0.3. This corresponds to a mean column density 
of hydrogen nuclei of (Nb) ~ 2-6 x 10 20 cm -2 , which com- 
pares well with the finding by the Copernicus satellite that 
there is a sudden jump in the H2 abundance along galac- 
tic sightlines wit h log A^h > 20.7, co nsistent with the onset 
of self-shielding |Savage et al.l fl977l ). It also agrees within 
a factor of a few with the critical extinction ^4y, C rit — 0.5 
derived by iKrumholz. McKee fc Tumlinsonl (|2008h for ide- 
alized spherical clouds. 
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Figure 5. (a) Evolution of the mass-weighted mean H2 abun- 
dance with time in runs nlOOO-ZOl (solid line), n300-Z01 (dashed 
line), nlOO-ZOl (dash-dotted line) and n30-Z01 (dotted line), (b) 
Evolution of the mass-weighted mean CO abundance with time 
in runs nlOOO-ZOl (solid line), n300-Z01 (dashed line), and nlOO- 
Z01 (dash-dotted line). In run n30-Z01, the mean CO abundance 
remains smaller than 10 — 11 for the whole period plotted. Reduc- 
ing the metallicity clearly has a much more pronounced effect on 
the CO than on the H2. 

3.4 Estimating the CO-to-H2 conversion factor 

It is of particular interest to understand how the distinct 
chemical changes we see as we decrease (^4v) affect the CO- 
to-H.2 conversion factor, Xco- The decrease of (ico)m by 
almost four orders of magnitude between (Ay) = 10 and 
(Av) = 1, is suggestive of a strong correlation between Xco 
and (Av), but does not by itself allow us to quantify the re- 
lationship between these two quantities. The sharp decrease 
in (ico)m as we decrease (Av) is reflected in a sharp de- 
crease in the column density of CO along most sightlines, as 
illustrated in Figure [S] This also suggests a strong correla- 
tion between Xco and (^4v), but still does not allow us to 
make a quantitative statement. 

To quantify the relationship between Xco and (Av), 
we need to be able to at least approximately compute 

where Wco is the velocity-integrated intensity of the J = 




Time (Myr) 

Figure 6. (a) Evolution of the mass-weighted mean H2 abun- 
dance with time in runs nlOOO-UVO (solid line), n300-UV0 
(dashed line), nlOO-UVO (dash-dotted line) and n30-UV0 (dotted 
line), (b) As (a), but for the mass- weighted mean CO abundance. 



1 — > line of 1 CO. Since we can easily measure Nh 2 along 
any desired line of sight through our simulations, the diffi- 
culty comes from the calculation of the integrated intensity. 
If the gas were optically thin in the 1 — > line, then it 
would be trivial to compute Wco, since the CO emissivity 
at any point would depend only on the local values of the 
temperature, density and chemical abundances, and would 
be effectively decoupled from the conditions elsewhere in the 
cloud. 

Unfortunately, in many of our simulations, some or all 
of the lines of sight are optically thick in the 1 — > line, 
making it far more challenging to model the CO emission 
accurately. When the gas is optically thick, it becomes nec- 
essary to account for line absorption, which can be ignored in 
the optically thin case. Furthermore, the opacity and emis- 
sivity in the line at any one point in the cloud can couple 
to the opacity and emissivity at other points in the cloud. 
Therefore, to model the CO emission entirely accurately, it 
becomes necessary to calculate the opacities and emissivities 
for all points in the cloud simultaneously. This is typically 
accomplished through the use of an iterative method (see 
e.g. the review of techniques in Ivan Zadelhoff et~alll 19921 ). If 
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Figure 7. (a) Evolution of the mass-weighted mean H2 abun- 
dance with time in runs nlOO (black solid line), n30 (black dashed 
line), nlOO-ZOl (red solid line), and n30-Z01 (red dashed line), (b) 
As (a), but for the mass- weighted mean CO abundance. Note that 
in run n30-Z01, the mean CO abundance remains smaller than 
10~ n for almost the whole of the period plotted, becoming barely 
visible in the plot only after t ~ 22 Myr. 



the geometry of the gas is relatively simple (e.g. slab sym- 
metry or spherical symmetry), then the resulting problem 
may still be solved without too much effort. In the present 
case, however, the density, temperature, chemical composi- 
tion and velocity structure of the gas are all highly inhomo- 
geneous, and solution of the coupled equations is a complex 
and computationally intensive problem that lies outside of 
the scope of this initial study. 

If we are prepared to tolerate a certain amount of un- 
certainty, however, then it becomes possible to estimate the 
CO-to-H2 conversion factor relatively easily, even in opti- 
cally thick gas. To compute an estimate for the conversion 
factor, which we denote as Xco.ost to distinguish it from 
the true conversion factor Xco, we use the following pro- 
cedure. We begin by selecting a set of independent lines of 
sight through our simulation results. We choose lines of sight 
that run parallel to the X-axis, one per resolution element. 
This choice is arbitrary; a different orientation would yield 
similar results. We then compute H2 and CO column densi- 
ties along each of these lines of sight. We next convert each 
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Figure 8. (a) CO column density in run nlOO at time t = t cn d, 
viewed along a line of sight parallel to the rr-axis of the simula- 
tion, (b) As (a), but for run n30. The contrast with run nlOO is 
striking, and illustrates how a relatively small change in the mean 
extinction can produce a large change in the appearance of the 
cloud. 
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of the CO column densities into an optical depth, no, which 
represents an estimate of the optical depth of the gas in the 
CO J — 1 — > transition. To perform this conversion, we 
make three major simplifying assumptions. We assume that 
all of the CO is in local thermodynamic equilibrium (LTE), 
and that the gas is isothermal, with a temperature equal to 
the weighted mean temperature of the gas, where we use the 
CO number density as an appropriate weighting function 



I], , k T(i,j,k)n C o(i,j,k) 



(11) 



J- mean — ^ — ^ / - - 7 \ 

We also assume that the CO linewidth Av is uniform, and is 
given by At; = 3kms -1 . Our choice for Av is motivated by 
the fact that our simulations have an RMS turbulent velocity 
of 5kms _1 , and that we would expect the one-dimensional 
velocity dispersion to be a factor of \/3 smaller than this, 
i.e. ctid = 5/\/3 ~ Skms" 1 . Given these as sumptions, th e 
relationship between Xco and no becomes (| Tielensl 1 20051 ^ 



Aioc gi 

Tw = o — 3 J° 

871-14) 9o 



1 — exp 



Nco 
Av 



(12) 



where Aio is the spontaneous radiative transition rate for the 
J = 1 — > transition, vio is the frequency of the transition, 
E\o = hv\a is the corresponding energy, go and gi are the 
statistical weights of the J = and J = 1 levels, respectively, 
and fo is the fractional level population of the J = level. 
We take values for Aio and vio from the LAMDA database 
(jSchoier et al.ll2005l ). In the case that T mC an = 10K, a typical 
temperature in the highly molecular gas in our simulations, 
this procedure yields no — 5 x 10 _17 Xco, i-e. in this case, 
a CO column density of approximately 2 x 10 16 cm -2 corre- 
sponds to an optical depth no = 1. 

To convert from no to Wco, we use the same pro- 
cedure as in the curve of growt h analysis presented in 
iPineda. Caselli fc Goodman! (j2008T ). We write the integrated 
intensity as 



Wco = 7b Aw 



20(f)df, 



(13) 



where lb is the brightness temperature of the line, Av is 
the linewidth, and /3 is the photon escape probability. We 
approximate /3 by assumin g that it is t he same as for a 
plane-parallel, uniform slab (ITielenal2005l ): 



/3(r) 



[1 - exp(-2.34-r)]/4.68T t^7, 



4t [in (t/VtF)] 



1/2 



r > 7. 



(14) 



Finally, since we have previously assumed that the CO is in 
LTE, we simply set T h = T mcan . 

Using Equations [13] and 1141 we can compute the inte- 
grated intensity for each of our lines of sight. Since we al- 
ready know the H2 column density along each line of sight, 
we could in principle associate a CO-to-tb conversion factor 
with each one. However, as the size of each individual zone 
in one of our simulations is much smaller than the spatial 
resolution of observations of all but the few closest GMCs, 
we feel it leads to a more meaningful value for Xco, est if we 
average over the whole cloud, since this is essentially what 
the limited resolution of the observations forces us to do 
obser vat ionally. 

We therefore average over all lines of sight to compute a 
mean intensity (W co ) for the datacube as a whole. Similarly, 
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Figure 9. Estimate of the CO-to-H2 conversion factor Xco, est, 
plotted as a function of the mean visual extinction of the gas, 
(Av). The simplifications made in our modelling mean that each 
value of Xqq est is uncertain by at least a factor of two. At 
(Ay) > 3, the values we find are consistent with the value of 
Xco = 2 X 10 20 cm~ 2 K —1 km -1 s dete rmined observationally 
for the Milky Way by iDame et al.l j200ll') . indicated in the plot 
by the horizontal dashed line. At (Ay) < 3, we find evidence for 
a strong dependence of Xco, est on (Ay). The empirical fit given 
by Equation [16] is indicated as the dotted line in the Figure, and 
demonstrates that at low (Ay), the CO-to-H2 conversion factor 
increases roughly as Xco, est Ay 3 ,J . It should also be noted 
that at any particular (Ay), the dependence of Xqq est on metal- 
licity is relatively small. Previous claims of a strong metallicity 
dependence likely reflect the fact that there is a strong depen- 
dence on the mean extinction, which varies as (Ay) tx Z given 
fixed mean cloud density and cloud size. 



we compute a mean H2 column density for the datacube, 
(Nh 2 ), by averaging over all lines of sight. Our estimate of 
the CO-to-H2 conversion factor is then simply the ratio of 
these two quantities: 



Xco,. 



(Nh 2 ) 

(w co )- 



(15) 



Using this procedure, we have computed this estimated 
CO-to-H2 conversion factor for each of our simulations, and 
have plotted the resulting values against (Ay) in Figure [9]- 
1101 In Figure [9] we show the results from all of our runs 
with non-zero UV backgrounds. In Figure 1101 we compare 
our results from the runs no UV with those from the corre- 
sponding runs with UV. The canonical value of Xco for the 
Milky Way, X CO = 2 x 10 20 cm~ 2 KT 1 km" 1 s (|Dame et all 
2001!) is indicated in both Figures by a horizontal dashed 
line. 

The first point to note from Figure [9] is that for our 
solar metallicity runs with (Ay) > 3, or in other words, for 
simulations with parameters that are a reasonable match 
to the properties of local GMCs, we find remarkably good 
agreement between our estimated conversion factor and the 
measured value. This level of agreement is probably fortu- 
itous to some degree, given the uncertainty in the obser- 
vational value and the larger uncertainties in our estimated 
values. Nevertheless, it does indicate that despite our simpli- 
fications, our estimates do seem to be successfully capturing 
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Figure 10. Estimate of the CO-to-H2 conversion factor -Xco, est; 
plotted as a function of the mean visual extinction of the gas, 
(Ay). In this Figure, we compare the values obtained for our solar 
metallicity runs with and without the presence of an ultraviolet 
background. At high (Ay), the presence of a UV background has 
little effect on our derived value of -Xco, est, but at low (Ay), 
CO photodissociation plays an important role in producing the 
observed Milky W ay value of -Xco = 2 X 10 20 cm~ 2 K _1 km -1 s 
l|Dame et al]|200ll) . indicated here by the horizontal dashed line. 



the basic physics and are probably trustworthy to within a 
factor of a few. 

A second important point to note from Figure [9] is that 
for any given (Ay), the CO-to-H.2 conversion factor shows 
little dependence on metallicity - changes to Z of an order 
of magnitude or more alter Xco.est by at most a factor of 
a few. It is therefore not true to say that Xco has a strong 
dependence on metallicity. Instead, Figure [9] demonstrates 
that the true dependence is on the mean extinction (Ay). 
For (Ay) > 3, there is no clear correlation between Xco, est 
and (Ay), but at (Ay) < 3, we see a strong correlation. 
The dependence of Xco, est on (Ay) can be described with 
the following, wholly empirical fitting function, indicated in 
Figure [9] by the dotted line: 



Xco,e 



2.0 x 10 20 

2.0 x 10 20 (Av/3.5)" 3 ' 5 



Ay > 3.5 
Ay < 3.5 



(16) 



where Xco, est has its standard units of cm -2 K" 1 km -1 s. 
This empirical fitting function is accurate to within a factor 
of a few over the whole range of (Ay) and Xco, est considered 
in the Figure. 

Of course, if we consider clouds with a fixed density and 
size, and merely vary the metallicity, then the strong depen- 
dence of Xco, est on the mean extinction implies a strong 
metallicity dependence, since (Ay) is proportional to metal- 
licity. On the other hand, if we vary the metallicity but also 
allow the mean density or the size of the cloud to vary to 
compensate, then we will find little or no metallicity depen- 
dence. 

FigurefTOldemonstrates the importance of UV photodis- 
sociation for producing the Galactic CO-to-tb conversion 
factor. In the absence of UV photodissociation, we get com- 
parable values of Xco, est for our (Ay) ~ 10 and (Ay) ~ 30 
runs, as one would expect, given that most of the gas in 




Mean visual extinction 

Figure 11. Velocity-integrated intensity of the J = 1 -> tran- 
sition in 12 CO, plotted as a function of the mean extinction of 
the gas. The symbols are the same as in our previous plots. Noise 
levels for extragalactic CO observations are typically of the order 
of 1 Kkms" 1 , so only the runs in the top right-hand quadrant of 
the plot produce enough CO to be readily detectable. 



these runs is well shielded from whatever UV is present, but 
we find significantly smaller values for Xco, est at lower mean 
extinctions. 

It is also of interest to plot the mean values of Wco for 
the various runs, as this gives a good indication of which of 
our simulations would correspond to CO-bright clouds, from 
which it would be easy to detect CO emission, and which 
would produce CO-faint clouds that would be difficult or 
impossible to detect in CO. In Figure ITT1 we therefore plot 
(Wco) against (Ay) for each of our simulations. Given that 
the noise level for observations of CO in Local Group galax- 
ies is typically of the order of 1 K km s _1 (e.g. iRosolowskvl 
120071 ) . it is clear that only the runs with (Ay) > 1 will pro- 
duce CO-bright clouds; runs with lower (Ay) will produce 
clouds that would be extremely difficult to detect with CO 
observations. Moreover, of the runs with (Ay) ~ 1, only 
runs nl000-Z003 and n300-Z01 produce detectable CO: the 
higher metallicity, lower density runs with (Ay) ~ 1 yield 
CO integrated intensities that would be undetectable. 

Finally, it is interesting to speculate how our results 
for Xco would be affected if we were to consider stronger 
UV field strengths. In clouds with a high mean extinction, 
we would not expect the increased UV to have much effect. 
In these runs, most of the CO is situated in gas that is 
already well shielded from the UV background, and even 
relatively large changes in the strength of the background 
will not lead to CO photodissociation becoming important 
in these regions. On the other hand, in runs with low mean 
extinction, that have most of their CO in gas that is not 
particularly well shielded from the UV background, a change 
in the strength of the background will have a much stronger 
effect. We would therefore expect the values we obtain for 
Xco in the (Ay) > 3 runs to have little sensitivity to the 
UV background field strength, and the values we obtain in 
the (Ay) < 3 runs to increase significantly with increasing 
UV field strength, leading to a steepening of the relationship 
between Xco and (Ay). We plan to test this prediction in 
future work. 
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4 DISCUSSION 

We can draw a number of immediate conclusions from the 
results that we have presented in the previous section. First, 
it seems plain that when we talk about the formation of a 
molecular cloud, we should draw a distinction between the 
formation of a cloud in which the hydrogen is primarily in 
molecular form and the formation of a cloud from which CO 
emission can readily be detected, since the conditions for the 
former are not the same as those for the latter. 

To form an H2-dominated cloud, the cloud must have 
a mean visual extinction of a few tenths of a magnitude 
(cf. iKrumholz. McKee fc Tumlinsonl lioosl . who find a crit- 
ical value of approximately 0.5 for a spherical molecular 
cloud). More importantly, however, the cloud must survive 
for long enough to convert most of its atomic hydrogen into 
H2. The time required is a strong function of noZ, the prod- 
uct of the mean density of the cloud and its metallicity, al- 
though the dependence is not linear, and the time required 
also depends upon the stren gth and nature of the turbulence 
|Glover fc Mac Lowll2007bl ; Milosavljevic et al, in prep.) 

On the other hand, to make a CO-bright cloud, we need 
a significantly higher mean extinction, (Ay) ^ 3, but do not 
require so much time, since the CO abundance in most of 
the gas comes into equilibrium within 1-2 Myr. In practice, 
given plausible cloud sizes, clouds that have a high enough 
mean extinction to be CO-bright will also be dense enough 
to form H2 relatively quickly, but the converse is not neces- 
sarily true. 

Second, the existence of a visual extinction threshold 
above which clouds become CO-bright provides a simple ex- 
planation for the observation that in low metallicity systems 
cloud masses derived from CO observations are significantly 
smaller than those derived from techniques th at do not de- 
pend on CO, such as infrared emission (e.g. Ilsraell 1 19971 ; 
iLerov et ai1l2007l . [20091 '). In these systems, the CO largely 
traces the regions of a cloud (or cloud complex) that have 
mean extinctions greater than 2, but does not trace the H2 in 
the surrounding envelope, which may be considerably more 
extensive and may contain a large fraction of the total mass 
of the cloud. 

Third, the fact that the CO fraction typically reaches 
equilibrium on a timescale comparable to or shorter than 
the crossing time of the cloud, and significantly shorter 
than the time required to assemble the cloud from warm 
atomic gas implies that CO emission will "switch-on" 
ra pidly during the assembly of the cloud, a s hypo thesized 
by lHartmann. Ballesteros-Paredes fc Berginl |200ll 'l. In the 
time it takes compressive flows or gravity to double or triple 
the mean extinction of a cloud with (Ay) ~ 1 its CO con- 
tent and luminosity can increase more than a hundredfold. 
The cloud will therefore quickly move from being effectively 
unobservable in CO to being readily observable. 

On the other hand, the long H2 formation timescale 
implies that molecular clouds need not have equilib- 
rium H2 abundances, even once they become CO- 
bright, contrary to what is often assumed (see e.g. 
Krum holz. McKee fc Tumlinsonl 120081 ; iMcKee fc Krumholzl 
How far from equilibrium the H2 fractions are will 
depend on the cloud properties (mean density, metallicity, 
rms turbulent velocity, etc.), the age of the cloud, and also 
on the range of densities considered within the cloud, since 



overdense regions will reach equilibrium much faster than 
underdense regions. 

Fourth, we would expect the molecular clouds (or re- 
gions thereof) in low metallicity systems that we observe to 
be CO-bright to also be systematically larger and/or denser 
than their counterparts in the Milky Way. This conclusion 
follows from the fact that if (Ay) > 3 is required for CO to 
form, and we decrease Z, then we must increase either no or 
L to compensate, since (Ay) oc noLZ. This conclusion also 
lends itself to a relatively simple observational test. If we 
observe a giant molecular cloud of size L in a low metallic- 
ity system (where L is the total size of the cloud, not of the 
CO-bright region), and compare it with a cloud of a similar 
size in the Milky Way, then the former should be system- 
atically denser than the latter. This systematic difference 
in density may be detectable by examination of the CO (2- 
1)/CO(1-0) or CO(3-2)/CO(l-0) line ratios, or the ratio of 
of HCN emission to CO emission, all of which will be larger 
in denser systems. 

Finally, our results also suggest a relatively simple ex- 
planation for the fact that GMCs observed in the Milky 
Way span only a small range in column densities (Blitz et al. 
120071 ). The requirement that (Ay) >3 for a CO-bright cloud 
implies that any cloud identified as a CMC must have a 
minimum mean column density of around 60 Mq pc -2 , since 
clouds with lower column densities will typically not contain 
much CO. 

The distribution of the mean column densities of GMCs 
above this minimum value depends upon the distribu- 
tions of the size and mean density of the clouds. How- 
ever, the cloud size distribution is typically found to be 
a shar ply decreasing function o f cloud size L. For ex- 
ample, iHever. Carpenter fc Snelil (|200lf ) find n(L)dL oc 
L" 3 2±01 dL for GMCs in the outer Galaxy. If the mean 
density of the GMCs spans only a small range of values, 
then the distribution of column densities n(E)dE will be a 
sharply decreasing function of E. If we combine this fact with 
the existence of a minimum column density threshold for a 
CO-bright GMC, then it implies that most observed GMCs 
will have column densities that are close to the threshold. 

This argument relies on the mean cloud density not 
being a strong function of the size of the cloud. It is pos- 
sibl e to show that this is a good a ssumption for at least 
the lHever. Carpenter fc S ncll (2001) cloud sample. Given a 
cloud size distribution n(L) oc ^- 3 - 2 ± 01 ] anc i a constant 
cloud mean density, then the cloud mass distribution is sim- 
ply 

n(M)dM = «(L)^dM, 

oc L- 3 ^^ 2 
oc M- (5 ' 2±0 - 1)/3 

oc M - 1Mm , (17) 

where we have assumed only that M oc L 3 , consis- 
tent with our adoption of a constant cloud mean density. 
For comparison, the cloud ma ss distribution inferred by 
IHever. Carpenter fc Snelil (|200lf ) for this sample of clouds 
was n(M) oc Af — 1 " 8±0 03 , demonstrating that even if the 
mean density of the clouds is not constant, it can at most 
be a very weak function of cloud size L. 

If this argument is correct, it suggests that most Galac- 
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l pc 

This value is a factor of th ree smaller tha n the value of 
170 Mq pc -2 determined by ISolomon et al . (1987) that is 
often taken to be canonical, but is in goo d agre ement with 
the value recently der ived by Hever et al.l (|2009f ) from their 
re-examination of the ISolomon et al.l clouds using new data 
from the Boston University-FCRAO Galactic Ring Survey. 

We note in conclusion that this argument makes no as- 
sumptions about the dynamical equilibrium of the clouds. 
This casts doubt on the frequently-made assumption that 
the small scatter in observed GMC column densities implies 
that the GMCs are in virial equilibrium. 



5 CAVEATS 

There are several caveats that should be borne in mind when 
considering our results. In our simulations, we use a simpli- 
fied chemical network for modelling the formation of H2 and 
CO, and adopt an approximate treatment of the effects of 
the ultraviolet background radiation field. When modelling 
the CO integrated intensity, Wco, we make additional ap- 
proximations, such as the assumption that all of the CO is 
in LTE, or that the CO linewidth has a fixed value. Finally, 
our simulations neglect the effects of self-gravity, which is 
clearly important in real molecular clouds. 

However, we are confident that these issues will not sig- 
nificantly affect our major results. In Paper I, we showed 
that our simplified chemical model produces results that 
agree exce llently with the standard UMIST astrochemi- 
cal model \Le Teuff, Millar, fc Markwickllioooh over a wide 
range of parameter space, provided that H, He, C and O are 
the only elements considered. Including other elements in 
the chemical model, such as sulphur, will worsen the agree- 
ment slightly, particularly in regions dominated by C and 
C + , but is unlikely to result in large changes in either the 
mean abundance of CO or the CO integrated intensity. 

Our simplified treatment of the effects of UV photodis- 
sociation will also introduce some error into our predicted 
CO abundances. Because the dense gas in a molecular cloud 
has a relatively small volume filling factor, we can be rea- 
sonably confident that our six-ray treatment of the radiation 
field will tend to overestimate its strength, as it will be more 
likely to miss shielding coming from dense regions that do 
not lie on lines of sight parallel to the coordinate axes, than 
to overestimate the effectiveness of the shielding provided by 
dense regions that do lie along these lines of sight. However, 
the main contribution to the photodissociation rate in any 
particular zone will come from the lines of sight with the 
lowest column densities. Therefore, if the lines of sight that 
run parallel to the coordinate axes have low column densi- 
ties, the fact that we miss the shielding from a few dense 
regions that do not lie along these lines of sight will not 
substantially alter the predicted photodissociation rates. 

The approximations that we have made in order to de- 
rive Wco are another potential source of error. Observation- 
ally, we know that not all of t he CO in molecular clouds is 
in LTE IjGoldsmith et alj|2008f ). and so our assumption that 
it is will cause us to overestimate the integrated intensity. 
Our assumption that all of the CO has the same tempera- 
ture is another source of error, although we have attempted 
to mitigate this by adopting a suitably-weighted mean tem- 



perature for the gas that gives rise to the CO emission. Our 
adoption of a fixed linewidth for the CO is also unrealistic: 
in reality, CO in lower metallicity or lower density systems 
will be distributed more intermittently, and so can be ex- 
pected to have a lower velocity dispersion than the gas as a 
whole. However, none of these sources of error will result in 
a large change in Wco, and moreover they will offset each 
other to some degree. 

Finally, our neglect of the effects of self-gravity is im- 
portant if one is concerned with the structur e of the cloud 
on the scales of individu al dense clumps (e.g. lKlessenll200ll ; 
iKainulainen et all 120091 ). However, the quantities that we 
have been dealing with in this paper are the results of averag- 
ing over the whole cloud, and changes in the cloud structure 
on scales of a few tenths of a parsec are unlikely to signifi- 
cantly affect these cloud-averaged values. Nevertheless, this 
is an issue that we intend to return to in the future. 



6 SUMMARY 

We finish our discussion here by summarizing the key re- 
sults of this paper. We have shown that the H2 abundance 
in turbulent molecular clouds is controlled primarily by the 
time taken to form the H2, and that it is relatively insensi- 
tive to the effects of UV photodissociation. Photodissocia- 
tion becomes significant for determining the H2 abundance 
only below a visual extinction threshold of a few tenths of a 
magnitude, in good agreement with previous observational 
and theoretical work. On the other hand, CO forms rapidly, 
but is strongly affected by photodissociation. The mean CO 
abundance falls off rapidly with decreasing mean extinction, 
particularly in clouds with (Ay) < 3. 

We have demonstrated that the CO-to-H2 conversion 
factor is also determined primarily by the mean extinction 
of the cloud, and that it is almost constant for (Ay) > 3, but 
falls off as Xco oc Ay 3 ' 5 for (Ay) < 3. Furthermore, only 
the clouds with visual extinctions greater than a few have 
sufficiently large CO integrated intensities to be detected 
by current observations. This therefore suggests a simple 
explanation for the discrepancy observed in low metallic- 
ity systems between cloud masses determined by CO obser- 
vations and those determined by non-CO tracers, such as 
IR emission. CO observations inevitably select CO-bright 
clouds that have mean extinctions large enough to place 
them in the regime where -Xco is approximately constant 
and equal to the Galactic value. On the other hand, obser- 
vations that do not rely on CO have no such bias and so 
select clouds that lie on the power-law portion of the rela- 
tionship. Therefore, CO observations find an Xco that does 
not vary significantly with metallicity, while other observa- 
tions find a strong dependence on metallicity. 

Finally, we have shown that if we combine the require- 
ment that (Ay) > 3 for a CO-bright cloud with the obser- 
vational fact that the GMC size distribution is a steeply 
decreasing function of size, then we are lead quite naturally 
to the prediction that all GMCs should have near-constant 
column densities. Therefore, contrary to what is often as- 
sumed, the small scatter in observed GMC column densities 
does not necessarily imply that the GMCs are in virial equi- 
librium. 
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